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Abstract — We present examples of agent-based and stochas- 
tic models of competition and business processes in economics 
and finance. We start from as simple as possible models, 
which have microscopic, agent-based, versions and macroscopic 
treatment in behavior. Microscopic and macroscopic versions of 
herding model proposed by Kirman and Bass diffusion of new 
products are considered in this contribution as two basic ideas. 
Further we demonstrate that general herding behavior can be 
considered as a background of nonlinear stochastic model of 
financial fluctuations. 

Keywords -agent based modehng; stochastic modeling; busi- 
ness models; financial market models. 

I. Introduction 

Statistically reasonable models of social systems, first 
of all stochastic and agent based, are of great interest 
for wide community of interdisciplinary researchers dealing 
with diversity of complex systems ||T|. Computer modeling 
serves as a technique in the quest for the relation between 
microscopic interactions of individual agents and macro 
dynamics of the whole system. Nevertheless, some general 
theories or methods that are well developed in the natural 
and physical sciences can be helpful in the development of 
consistent micro and macro modeling of complex systems 

Our approach in the modeling of the financial markets, 
using the nonlinear stochastic differential equations, was 
primarily based on the empirical analysis of financial market 
data. The proposed equations have power law statistics very 
similar to the ones observed in the empirical data (the broad 
review of our works in this field is given in |3|). Though 
all of this work was done by relying on the macroscopic 
phenomenological reasoning, the microscopic, interaction 
based, grounding is of the utmost interest to both scientific 
and policy-making community |4|, fS), ID. The development 
of the macroscopic treatments for the well established agent 
based models would be the most consistent approach in the 
attempt to understand the essential correspondence of micro 
and macro behavior. For this task one should select the 
simple agent based models, which would have an expected 
macroscopic description. In this contribution we present a 



few examples from the agent based modeling in the business 
and finance while showing that the example have useful and 
informative macroscopic treatments. 

Kirman's ant colony model [7J is an agent-based model 
used to explain the importance of herding inside the ant 
colonies and economic systems (see the later works by 
Kirman and other authors, which develop on this idea 
Q, im, ||9|)- The analogy can be drawn as human crowd 
behavior is ideologically and statistically similar in many 
senses. On our website, ifTOl . we have presented interactive 
realizations of the original Kirman's agent-based model (see 
ifTTI ') and of its stochastic treatment by Alfarano et al. (SI 
(see ifTl). 

The diffusion of new products is one of the key problems 
in marketing research. The Bass Diffusion model is a very 
prominent model related to this problem. This model is 
formulated as an ordinary differential equation, which might 
be used to forecast the number of adopters of the new 
successful product or service |13|. There were suggestions 
that such basic macroscopic description in marketing re- 
search can be studied using the agent-based modeling as 
well |[T4|. Thus it is a great opportunity to explore the 
correspondence between the micro and macro descriptions 
looking for the conditions under which both approaches 
converge. The Bass Diffusion model is of great interest for 
us as representing very practical and widely accepted area of 
business modeling. Web based interactive models, presented 
on the site [15] serve as an additional research instrument 
available for very wide community. 

Further we follow the works by Alfarano et al. fS^, f9\ and 
introduce our own model modifications in order to obtain 
more sufficient agent-based models of financial markets, 
which would have an alternative macroscopic description 
in the terms of Stochastic Calculus. We also simplify the 
Kirman model to the unidirectional case and show that it 
can be seen as microscopic analog to the Bass diffusion 
model (the interactive applet of this model is also available 
on our website HH). 

In Section II we briefly introduce our online publishing 
and collaboration platform, while the further Sections are 



devoted: III to the macroscopic treatment of Kirman's model, 
IV to the agent base modeling of Bass diffusion, V to 
the modeling of financial markets by nonlinear stochastic 
differential equations, VI to the application of Kirman's 
model to the financial markets, VII to the multifractal 
behavior of return series, VIII to the statistics of bursts 
generated by nonlinear SDE and finally we present some 
conclusions and discussion on future work. 

II. Web platform 

Our web site flU] was setup using WordPress webloging 
software. The setup pays to be user-friendly, powerful and 
easily extensible web publishing platform, which with some 
effort can be adapted to the scientist's needs. There is a 
wide choice of plugins, which enable convenient usage of 
equations (mostly using LaTeX). While during the setup we 
found that bibliography management plugins were lacking 
at the time. 

To accommodate our needs for the equations we have 
worked on improving WP-Latex plugin (available from 
|17|). Namely we have introduced a possibility to write 
equations in both inline and ordinary math modes. Imple- 
mented the equation labeling, numbering and referencing. 
And finally fixed some noticeable problems with vertical 
placement of the inline mode equations. 

Another important task was to implement bibliography 
management and citations. For this cause we have used 
the bibtexParse PHP code (available from ifTSi ) to setup 
BiBTeX backend. From this point on we have written our 
own original PHP code to link between bibtexParse, our 
database and WordPress. By using this plugin we can now 
easily manage and present our own papers (ex. generate 
our own bibliographies), papers we have read (tag them 
with keywords, write our own comments and etc.) and also 
communicate with the visitors using numerous citations. 

Interactive models themselves are independent from the 
publishing framework. Most of them were implemented 
using prominent Java applet technology. Some of the ap- 
plets were created using multi-paradigm simulation software 
AnyLogic |19], while the others were programmed from 
scratch using Java programing language 1201 . AnyLogic was 
used in the most of agent based scenarios as it is a very 
convenient tool for agent based modeling, while programing 
from scratch gave us more control over the applets behavior 
needed while doing stochastic modeling. 

Either way by compiling appropriate files one obtains Java 
applets, which can be included in to the articles written 
using WordPress. This way articles become interactive - 
visitor can both theoretically familiarize himself with the 
model and test if the claims made in the post describing 
model were true. This happens in the same browser window, 
thus the transition between theory and modeling appears to 
be seemless. Due to the fact that models are implemented 
as Java applets al of the numerical evaluation occurs on 



client machine, while the visitor must have Java Runtime 
Environment installed, and server load stays minimal. The 
requirement for JRE might appear to be cumbersome, but the 
technology is somewhat popular and freely available from 
Oracle Corp. 

One of the goals of developing these models on the web 
site ifTSl was to provide theoretical background of Bass 
Diffusion model and discuss practical steps on how such 
computer simulations can be created even with limited IT 
knowledge and further applied for varying purposes. Thus, 
we have targeted small and medium enterprises to encourage 
them to use modem computer simulation tools for business 
planning, sale forecasting and other purposes. 

Consequently computer models and their corresponding 
descriptions published at the | fT5| provide a relatively easy 
starting point to get acquainted with computer simulation 
in business. The published content enables site visitors 
to familiarize themselves with these models interactively, 
running the applets directly in a browser window, changing 
the parameter values and observing results. This significantly 
increases accessibility and dissemination of these simula- 
tions. 

III. Extended macroscopic treatment of 
Kirman's model 

There is an interesting phenomenon concerning behavior 
of ant colony. It appears that if there are two identical food 
sources nearby, or two identical paths to the same food 
source (the experiment done by Pasteels and Deneubourg 
IfZTI . 1221), ants exploit only one of them at a given time. 
Evidently the food source which will be used at a given 
time is not certain. It is so as switches between food sources 
occur, though the food sources, or paths, remain the same. 

One could assume that those different food sources are 
different trading strategies or, if putting it simply, the actions 
available to traders. Thus, one could argue that speculative 
bubbles and crashes in the financial markets are of similar 
nature as the exploitation of the food sources in ant colonies 
- as quality of stock and quality of food in the ideal case 
can be assumed to be constant. Thus, model |7 1 was created 
using ideas obtained from the ecological experiments BTI . 
1 22 1 can be applied towards the financial market modeling. 

Kirman, as an economist, actually developed this model 
as a general framework in context of economic modeling 
(see Q, ll23l and his other works). Recently his framework 
was also used by other authors who are concerned with 
the financial market modeling (see jS), ||9l). Thus basing 
ourselves on the main ideas of these authors and our pre- 
vious results in stochastic modeling (see |3|) we introduce 
specific modifications of Kirman's model providing a class 
of nonlinear stochastic differential equations |24| applicable 
for the financial variables. 

Kirman's one step transition probabilities might be ex- 



pressed in the following form 121, 

p{X ^ X + 1)={N - X) (ai + hX) At, (1) 
p{X ^ X - 1)=X {c72 + h[N - X]) At, (2) 

where X is a number of agents exploiting the chosen trading 
strategy (the one used to describe system state), while N 
is a total number of agents in the system (thus the other 
trading strategy is used by the N — X agents). In the above 
the original Kirman's approach was extended by introducing 
fixed event time scale At by replacing the original models 
individual decision Si (JiAt and herding [1 — 5) ^ hAt 
parameters. Later we will need a more general assumption 
that parameters cr and h may depend on X and N, but for 
now we omit it. 

Note that the transition probabilities ([T]i and Q describe a 
scenario where the interaction among agent groups depends 
on the overall number of agents in alternative state. Such a 
choice makes the transition rates non-extensive, the connec- 
tivity between agent groups increases with the number of 
agents N. The herding interactions have a global character. 
Opposite scenario - extensive one will be also used further 
in this paper. 

The lack of memory of the agents is the crucial assump- 
tion to formalize the population dynamics as a Markov 
process. Furthermore to describe the aforementioned dy- 
namics in a continuous time we will need to obtain the 
transition rates, transition probabilities per unit time, which 
for continuous x — X/N may be expressed as 



By introducing custom functions 



Tr+{xMl-x)Q + hx), 
7r-{x)=x(^^+h[l-x]). 



(3) 
(4) 



Here the large number of agents N is assumed to ensure 
the continuity of variable x, which expresses the fraction 
of agents using the selected trading strategy, X. Relation 
between the discrete transition probabilities, ([T]i and and 
continuous transition rates, ^ and (|4]i, should be evident; 

p{X X ±1) ^ N^TT^ {x)At. (5) 

One can compactly express the Master equation for the 
system state probability density function, U!{x, t), by using 
one step operators E and E^-*^ (see ll25l for a details on this 
formalism) as 



dtu}{x,t) =N^ {{^ - l)['K- {x)Lo{x,t)] + 
+ (E-i - l)[Ti+{x)oj{x,t)]] 



(6) 



By expanding E and E^^ using the Taylor expansion (up 
to the second term) we arrive at the approximation of the 
Master equation 



dtuj{x,t) ^-Nd^[{TT+{x) - TT-{x)}uj{x,t)\ + 
+ \dl[{^+{x)+^-{x)}uj{x,t)]. 



A{x)^N{t:+{x) - TT^{x)} = cri(l - x) - a2X, (8) 
D{x)^TT+{x) +TT^{x) = 2hx{l -x) + 

(9) 



I ^ I ^2 



one can make sure that the (j6]) is actually a Fokker-Planck 
equation; 

dMx,t) = -d,[Aix)u;ix,t)] + ^dl[D{x)uj{x,t)]. (10) 

Note that in the limit of large N one can neglect individual 
behavior terms in the D{x). The above Fokker-Planck 
equation was first derived in a slightly different manner in 
the ISJ. 

It is known (for details see f26l) that the Fokker-Planck 
equation can be rewritten as Langevin equation, or in other 
words stochastic differential equation, 

dx^A{x)dt + ^D{x)dW = 

= [cri(l -x)- (J2x]dt + ^2hx{l - x)dW, (11) 

here W stands for Wiener process. This step was also present 
in the fSl. 

In Fig. [T] we show that the statistical properties ob- 
tained from the agent based model, defined by transition 
probabilities ([T]) and Q, match statistical properties of 
the solutions of ( [TT] i. Thus the approximations done while 
deriving the Langevin equation for population fraction are 
valid. Interestingly enough we have obtained agreement with 
not so high number of agents - N — 100. 

Note that the method used to derive Eq. (Ill gives us an 
opportunity to consider parameters ai, G2, h dependent on 
the variable x and N. We will need this generalization in the 
further elaboration on various applications. From our point 
of view, the general form of SDH ( fTT) derived from the very 
basic agent based herding model provides a wide choice of 
opportunities in consistent micro and macro modeling of 
complex social systems. 

IV. Agent based model for the Bass Diffusion 

The Bass Diffusion model is a tool to forecast the diffu- 
sion rate of new products or technologies ifTSl . Mathemati- 
cally it is formulated as an ordinary differential equation 

h 



dtX{t)=[N - X{t)] 
X{0)^0. 



N 



X{t) 



(12) 
(13) 



(7) 



where X{t) denotes the number of consumers at time t, N 
can be seen as the market potential, being a starting number 
of the potential consumers (agents), a is the coefficient 
of innovation, the likelihood of an individual to adopt the 
product due to influence by the commercials or similar 
external sources, h is the coefficient of imitation, a measure 
of likelihood that an individual will adopt the product due to 
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Figure 1. Agreement between statistical properties of population fraction, x, (a) probability density function and (b) power spectral density, obtained 
from stochastic (red and blue curves) and agent-based (green and magenta squares) models. Two qualitatively different model phases are shown: red curve 
and green squares correspond to herding dominant model phase (cri = (T2 = 0.2, h = 5), while blue curve and magenta squares correspond to individual 
behavior dominant model phase (cri = cr2 = 16, h = 5). Agent based model results obtained with A'^ = 100. 



influence by other people who akeady adopted the product. 
This nonlinear differential equation serves as a macroscopic 
description of new product adoption by customers widely 
used in business planning |14|. 

The agent based approach to the same problem is related 
with modeling of product adoption by individual users, or 
agents. One can simulate diffusion process using computers, 
where individual decisions of adoption occur with specific 
adoption probability affected by the other individuals in the 
neighborhood. It is easy to show that Bass diffusion process 
is a unidirectional case of the Kimian's herding model Q. 
Indeed, let us define x{t) in the same way as in previous 
section x{t) — X{t)/N, then the potential users will adopt 
the product at the same rate as in Kirman's model agents 
switch from one state to another 



7r+(x) = (l 



X — 



(7 

N 



h 

— X 

N 



(14) 
(15) 



The form of ( fTSj ) should be self explanatory - in case of the 
product diffusion agent should not be allow to withdraw 
from the consumer state, thus this transition probability 
should be forced to equal zero. 

The mathematical form of ( [l4| l is not as evident, note 
that we have substituted h with (compare with the 
original model transition probability ([3]l), and needs further 
discussion. Mathematically this substitution can be backed 
by the need for the stochastic term to become negligible 
in the limit of large N. In the modeled market terms 
this substitution means an introduction of the interaction 
locality - namely it is an assumption that each individual 
communicates only with his local partners (epidemic case). 

One can compare the expression of the transition prob- 
ability, ( 14 1, with the adoption probabilities of the Linear 



and GLM models of Bass Diffusion discussed in lIZTl . The 
match in expressions is clear in the small time step limit. 
At -> 0. 



In case of the transition rates ( [T4] i and ([TSj the macro- 
scopic description functions, namely drift, A{x), and diffu- 
sion, D{x), become 

A{x)=NTr+{x) = {l-x){<7 + hx), (16) 

(17) 



D{x)=7T+{x) = -^^ ^ (a + hx) 



N 

In the large market potential limit, N ^ 1, D{x) becomes 
negligible and thus one can consider the obtained equation 
to be equivalent to the Bass Diffusion ordinary differential 
equation ( [T2| instead of the stochastic differential equation. 
This serves as a proof that Bass Diffusion is an unidirec- 
tional epidemic case of Kirman's herding model. Though 
this simple relation looks straightforward, we derive it and 
confirm by numerical simulations in fairly original way. 

In Figure |2] we demonstrate the correspondence between 
the Bass Diffusion model (macroscopic description) and 
unidirectional Kirman's herding model (microscopic de- 
scription). Both, agent based and continuous, descriptions 
of the product adoption, AX, converge while the market 
potential, N, or the selected observation time interval, r, 
become larger 

V. Nonlinear stochastic differential equation 

AS A MODEL OF THE FINANCIAL MARKETS 

Earlier we have introduced a class of non-linear stochastic 
differential equations (SDEs) providing time series with 
power-law statistics, and most notably reproducing 1//^ 
spectral density, ||28]| . Il29l . Il30l . The general form of the 
proposed class of Ito SDEs is 



dy= 77 - 



(18) 



here y is the stochastic process exhibiting power-law statis- 
tics, f] is the power-law exponent of the multiplicative 
noise, while A defines the exponent of power-law probability 
density function (PDF), and is a Wiener process (the 





Figure 2. Comparison of the product adoption per observation interval, AX/r versus t, obtained from the macroscopic description by the Bass Diffusion 
model, jl2^ , (red line) and the microscopic description using the unidirectional Kirman's model, jl4| , (blue points). The models tend to converge when 
time window, r, or market potential, N, become larger: (a) = 1000, r = 0.1; (b) N = 1000, r = 1; (c) Af = 10000, r = 0.1; (d) N = 10000, 
T = 1. Other model parameters were the same for all subfigures and were as follows a = 0.01, h = 0.275. 



Brownian motion). Note that SDE ([TSj is defined in the 
scaled time, tg = cr^t, where erf is the scaHng parameter. 
Empirically we have determined that af = 1/6 • lO^^s^^ is 
appropriate in terms of the return model proposed in lISTl . 



From the SDE ( 18 i follows that the stationary probability 



density function (PDF) of this stochastic process is power- 
law, po{y) ^ y^^, with the exponent A |26|. While in Refs. 
and later more precisely in ll24l it was shown that the 



time series obtained while solving SDE ( 18 1 have power-law 
spectral density 



1 



/? = 1 + 



A-3 
2(^7-1) 



(19) 



Note that exponent of spectral density, /3, is defined only for 
77 7^ 1. In case of = 1 the SDE ( fTSj l becomes identical to 
the geometric Brownian motion. 

Power law statistics of the signal y obtained by solving 



SDE ( 18 1 and exponents A, fi are defined for large y values. 
Thus one has to introduce the diffusion restriction terms 
in the limit of small y values when attempting to solve 



SDE (18 1 or applying it in a stochastic modeling. There 
is a wide choice of restriction mechanisms adjustable to the 
needs of real systems with negligible influence on the power 
law exponents. We have introduced a term of additive noise 



while attempting to model the absolute return ||3T1 

d2/= (^~^{l + y^y-^ydt, + {l^y^Y^AWs. (20) 

In such case the stationary probability density function of 
the SDE ([T8]l is a g-Gaussian (see Q, Il3l1 ) 



Pxiv) 



r(A/2) 



(21) 



y^r(A/2-i/2) Vl + y^ 

While modeling the trading activity |f33l we have used the 
exponential diffusion restriction for small values of variable 

y - ymin 



dy = 



/ 2 2 V r 



y 



^"^-hlt + y^dW. (22) 



Equation ( 22 1 has a very general form, which includes the 
well known models applicable to financial markets such 
as the Cox-Ingersoll-Ross (CIR) process or the Constant 
Elasticity of Variance (CEV) process ||34| 

dy = ^lydt + y''dVK, (23) 

where ^ — {rj ~ l)ymin ^ ^^^^ general cases of the 

SDE 

The class of equations based on SDE ( [T8] l gives only a 
general idea how to model power-law statistics of trading 



activity and return in the financial markets. The problem is 
to determine the parameter set A and 77 in a way giving the 
empirical values for the A and /3. The task becomes even 
more complicated if one considers the more sophisticated 
behavior of the spectral density - power spectral densities 
have not one, but two power-law regions with different 
values of f3. In the series of papers lISTl . Il33l . Il35l we 
have shown that trading activity and return can be modeled 



by a more sophisticated version of the SDE than (I81 now 
including the two powers of the noise multiplicativity. In the 
case of return instead of Eq. ( [20) l one should use 

2/ {eVT 



r 



1)^ 



(24) 



here e divides the area of diffusion into the two different 
noise multiplicativity regions to ensure the spectral density 
of |y| with two power law exponents. 

The proposed form of the SDE enables reproduction of 
the main statistical properties of the return observed in the 
financial markets. Similarly one can deal with a more sophis- 
ticated model for the trading activity lf33l . This provides an 
approach to the financial markets with behavior dependent 
on the level of activity and exhibiting two stages: calm 
and excited. Equation ( [24| models the stochastic return y 
with two power-law statistics, namely the probability density 
function and the power spectral density, reproducing the 
empirical power law exponents of the return in the financial 
markets. 

VI. KiRMAN'S MODEL AS A MICROSCOPIC APPROACH 
TO THE FINANCIAL MARKETS 

The drawback of the stochastic models is a lack of direct 
insights into the microscopic nature of replicated dynamics. 
Bridging between microscopic and macroscopic approaches 
is needed for better grounding of stochastic modeling. 

Top-down approach, namely starting from the stochastic 
modeling and moving towards the agent based models, 
seems to be a very formidable task, as the macro-behavior 
of complex system can not be understood as a simple su- 
perposition of varying micro-behaviors. While in the case of 
sophisticated agent based models 1361 bottom-up approach 
provides too many opportunities. But there is selection of 
rather simple agent based models (ex. |7|), whose stochastic 
treatment can be directly obtained from the microscopic 
description fE]. 

Here we consider an opportunity to generalize Kirman's 
ant colony model [7| with the intention to modify its 
microscopic approach to the financial market modeling fS\ 
reproducing the main stylized facts of this complex system. 
In the Section III we have already introduced Kirman's 
ant colony model, proposed its generalization and derived 
stochastic model for the two state population dynamics. 



As Kirman's model considers the two available agent 
states one must define two types of agents acting inside 
the market in order to relate Kirman's model to financial 
markets. Currently, the most common choice is assuming 
that agents can be either fundamentalists or noise traders 
136J. 

Fundamentalists are assumed to have the fundamental 
knowledge about the market, which is assumed to be quanti- 
fied by the so-called fundamental price, Pf{t), of the traded 
stock. By having this knowledge they can make long term 
forecasts on a notion that infinitely long under-evaluation or 
over-evaluation of the stock is impossible - the market in 
some point in the future will have to set a fair price on the 
stock. Thus their excess demand, which is shaped by their 
long term expectations, is given by |8| 



Df{t) = Nf{t)\n 



Pfjt) 
Pit)' 



(25) 



here Nf (t) is a number of fundamentalists inside the market 
and P{t) is a current market price. As long term investors 
fundamentalists assume that P{t) will converge towards 
Pf{t) at least in a long run. Therefore if Pf{t) > P{t), 
fundamentalists will expect that P{t) will grow in future 
and consequently they will buy the stock {Df{t) > 0). In 
the opposite case, Pf{t) < P{t), they will expect decrease of 
P{t) and for this reason they will sell the stock (Df{t) < 0). 

The other group, the noise traders are investors who 
attempt estimate the stocks future price based on its re- 
cent movements. As there is a wide selection of technical 
trading strategies, which are used to analyze stocks price 
movements, one can simply assume that the average noise 
traders demand is based on their mood, ^(i), |i8j 



(26) 



here Dc{t) is a total excess demand of noise trader group, 
Nc{t) is a number of noise traders inside the market and tq 
can be seen as a relative noise trader impact factor. 

Price and, later after a brief derivation, return can be 
introduced into the model by applying the Wakassian sce- 
nario. One can assume that trading in the market is cleared 
instantaneously to set a price, which would stabilize the 
market demand for a given moment. Thus the sum of all 
groups' excess demands should equal zero: 



Df{t)+D,,{t)=Nf{t)\n 



Pfjt) 
Pit) 



- roN,{t)at) 



=0, 

P(t)=P/(t)exp 



Ncit)^,.. 

"m)^^\ 



(27) 
(28) 



where without loosing generality one can assume that fun- 
damental price remains constant, — Pj- 



Consequently the return, which is defined as logarithmic 
change of price, in the selected time window T is given by: 



'•{t)=\nP{t)-\nP{t-T) = 



1 - x{t) 



x{t - T) 
l-x{t- T) 



at - T) 



,(29) 



where we have set that 



X and 



Nf{t) 



1 



according to the notation introduced in Section |lll] Alfarano 
et al. |8| simplified the above by assuming that x{t) is 
significantly slower process than obtaining adiabatic 
approximation of the return 



1 — x(t) 



(30) 



where = ^{t) — ^(t — T). If C(t) is modeled using spin- 
noise model, as in |8|, then the middle term. 



l-x{t) 



can be 



seen as an absolute return. 

Using Ito formula for variable substitution ll26ll in SDE 
we obtain nonhnear SDE for the y{t) = y^^^ 







dy={ai-y[a2-2h])il 



-y)dW. (31) 



Agreement between the agent based Kirman's model applied 
towards financial markets using the ideas discussed above 
and the new stochastic model for y, ( [ST) , is demonstrated 
in Fig. [3] 

Note once again that the actual derivation, and thus, the 
final outcome, does not change even if cti, (T2 or h are the 
functions of either x or y. Therefore, one can further study 
the possibilities of the obtained stochastic model, ( [3T] i, by 
checking different scenarios of cti, (T2 or h being functions 
of either x or y. Nevertheless, the most natural way is to 
introduce a custom function r(y) to adjust the inter-event 
time according to the system state. From the financial market 
point of view this can be seen as introduction of variability 
of trading activity based on the return. 

We have chosen the case when h and a2 are functions 
of y, namely we make the substitutions, a2 and 
h TMi Kirman's model transition probabilities, 

([T} and (|2]i, and stochastic model for y, pT) . To further 
simplify the model we can introduce scaled time, tg = ht, 
and make related model parameter transformations, £i — j^- 
By making these substitutions we arrive at 



dy 



si +y- 



£2 



T{y) 



{i + y)dts 



2y_ 
T{y) 



(1 + y)dWs 



(32) 

where Wg is appropriately scaled Wiener process. Note 
that we left cri, and consequently ei, independent of y 
on purpose as one could argue that individual behavior of 
fundamentalist trader should not depend on the observed 
returns as he is a long term investor uninterested in the 
momentary fluctuations of the market mood. 



Note that absolute return, y, defined in Eqs. ( 3 1 1 and ( 32 1 



serve as a measure of volatility in the financial markets. It is 
known that volatility has long-range memory and correlates 
with trading activity and has probability density function 
with power law tail [37l. We are particularly interested in 
the case of T{y) = y~". This selection is defined by the fact 
that trading activity has positive correlation with volatility 



and the class of SDE (18 1 is invariant regarding power-law 
variable transformation, see [301 . In such case the obtained 
stochastic differential equation, Eq. p2| , in the limit of y ^ 
1 is very similar to the stochastic models discussed in the 
Section [V] 

In the aforementioned Umit of ?/, y ^ 1, we can consider 



only the highest powers in Eq. (32 1. In such case Eq. (32i 
is reduce to the 



dy = (2 - £2)2;'+"dt, + v/2y3+"dW^,. 

The direct comparison of Eqs. ( fTSl and ( [33] l yields 

3 + a , 



(33) 



(34) 



Consequently we expect that the stochastic process y defined 



by Eq. (33 1 will have the power law stationary probability 
density function. 



-£2 — a— 1 



p{y) ^ y 

and also a power law spectral density, 

£2 + a -2 



1 



1 



a 



(35) 



(36) 



where we have used the relation between model parameters. 



Eq. (34 1 



While if we linearize drift function of Eq. ( 3 1 1 with the 



respect to the absolute return, y, namely set £2 = 2, we 
would obtain a stochastic differential equation (once again 
in the limit y ^ 1) 



dy = £iydts + v/2y3+"dW^, 



(37) 



similar to the generalized CEV process, which was consid- 
ered in Ii38i 

dy = aydt + by'^dW. (38) 



In I3OI the latter was noted to be a special case of Eq. 
(22 1 with exponential restriction of diffusion applied. The 
comparison with this special case is important on its own 
as this equation generalizes some stochastic models used in 
risk management. Theoretical prediction of PDF and spectral 
density for y defined by Eq. (|37]l, is given by L38i 

p{y) - y-'-'', (39) 

^ = 1 + ^, (40) 
where we have used the previously obtained relation between 



model parameters, Eq. (34i. 

In the Figure |4] we show that the theoretical predictions 
discussed in this section are valid and that they enable the 




Figure 3. Agreement between statistical properties of y, (a) probability density function and (b) power spectral density, obtained from the stochastic (red 
and blue curves) and agent-based (green and magenta squares) models. Two qualitatively different model phases are shown: red curve and green squares 
correspond to herding dominant model phase (a\ = (T2 = 0.2, h = 5), while blue curve and magenta squares correspond to individual behavior dominant 
model phase {a\ = a2 = 16, h = 5). Agent based model results obtained with N = 100. 



reproduction of different spectral densities and probability 
density functions. 

Note that while the stochastic model based on herding 
behavior of agents appears to be too crude to reproduce 
statistical properties of financial markets in such details 
as the stochastic model driven by the Eq. p4| , which is 
heavily based on the empirical research, it contains very 
important long range power law statistics of the absolute 
return. Obtained equations are very similar to some general 
stochastic models of the financial markets I.24J . [38 1 and 
thus, in future development might be able to serve as a 
microscopic justification for them and maybe for the more 



sophisticated model driven by the Eq. ( 24 1 



It is possible to extend agent based model by introducing 
additional agent groups or splitting old ones. Let us assume 
that chartist agents may disagree in their expectations and 
thus divide into pessimists and optimists. Therefore it is nat- 
ural to introduce three agent groups (see Fig. |5]l interacting 
among themselves. Our first attempts in this direction proves 
that in case of the three agent groups (as shown in Fig. |5]l, 
when the herding parameter hcc ^ hcf, might confirm the 
expectation of a more complex behavior exhibiting fractured 
power spectral density of absolute return. More detailed 
study of such approach in comparison with macroscopic 



modeling by SDE ( 24 1 is ongoing 



VII. MULTIFRACTAL BEHAVIOR OF RETURN SERIES 

In the last few decades it was noted that many natural 
phenomena have very complex intrinsic structure, which 
has a very specific scaling properties. This notion was 
generalized as fractal framework ll39l . Later it was also 
noted that the scaling properties of some processes exhibit 
even more complex scaling behavior - namely they appeared 
to have features of the multiple fractals. Few examples of 
such phenomena include geoelectrical processes ||40|, human 
heartbeat |41 J and gait |42). The financial market time series 
apparently are also of the multifractal nature [43 J , [441 . P31 . 



There are few established methods to detect multifractal 
time series and two very prominent methods. One of them 
is generalized height-height correlation function method 
(GHHCF) and multifractal detrended analysis method (ME- 
DEA). In our previous approaches B6ll . BtII we have used 
the GHHCE method, so let us in this contribution to rely on 
the ME-DEA method. 

To start with the multifractal analysis of the time series, 
Uk, we have to obtain the profile of the time series, Y^: 



(41) 



Next we have divide the Yk series into equally sized and non 
overlapping segments. Thus if our segments are of the size s, 
then we will have iVg = int(A^/s) segments (here N is the 
length the series, while int(. . . ) is a function which takes an 
integer part of the argument). Eor the most of the segment 
sizes some of the data will be lost, in order to account for 
it one might want to take another set of segments, but now 
splitting from the end of the series. 

Eurther, one has to determine the trends in the obtained 
segments. Generally this can be done using varying poly- 
nomial fits, but linear fits in the most cases are more than 
enough. After the trends, Y , are known one has to evaluate 
how well the trend fits the actual series: 

= - [^('^-l)s+^ - Mi)]". (42) 

1=1 

1 * 2 

X! [vN-(y-N,)s+i - yv{i)\ ■ (43) 



The Eq. ( [42] i holds for segments — l,...,Ns, while 
the Eq. (|43[) should applied towards segments ly = Ng + 
l,...,2Ns. Einally one has to average over all segments 
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Figure 4. Statistical properties, (a) and (c) - probability density function, (b) and (d) - spectral density, of the time series obtained by solving Eq. |32j 
(colored squares). Fits are provided by the theoretical predictions made in this section, (a) and (b) are fitted by using l |35( and |36j, (c) and (d) are fitted by 
using \39) and j40) , (curves of corresponding colors). Model parameters for the (a) and (b) were set as follows: a = 1, ei = 0, £2 = 0.5 (red squares), 1 
(green squares), 1.5 (blue squares) and 2 (magenta squares). Model parameters for the (c) and (d) were set as follows: ei = £2 = 2, a = (red squares), 
1 (green squares) and 2 (blue squares). 




Figure 5. The general case of the three groups of interacting agents: / - fundamentalists, c+ - chartists optimists, c_ - chartists pessimists, hij are 
herding terms, while a^, 6j and stand for individual transitions in the direction of the arrow. 



using 

^.(■^) = WE[^'(-)]^| ' (44) 

here q stands for generalized coefficient, which is the one 
enabling us to recover multifractal features it is also the only 
difference from the original detrended fluctuation analysis 
(DFA) method 148J. Note that in case of g = 2 the Fg{s) is 



the same as the one in the original DFA method. 

All that is left is to determine is the power law trend, h{q), 
of the Fq{s). These trends, h{q), are also frequently named 
the generalized Hurst exponents. If the Hurst exponents are 
different for different q, which can be any real number, then 
the signal can be seen as multifractal. In the opposite case 
or if the variation is negligible, time series can be assumed 
as monofractal. For more details on the MF-DFA method 



see 

In Figure |6] we show that the stochastic differential 
equations obtained for the modeHng of financial markets 
and derived from the Kirman's agent based model have 
broad multifractal spectra. The curves capture a region of 
the Brownian motion, h{q) w 1.5, and a region of long 
range memory, h{q) « 1. Note that in case of a = 1 and 
a = 2 (green and blue curves) h{2) = 1, which can be seen 
as a proof that the obtained time series posses the long term 
correlated (have so-called long range memory) behavior, 
while for a ~ interim behavior between the Brownian 
motion and long range memory is observed, 1 < h{2) < 1.5. 



VIII. 



Statistics of bursts generated by 
nonlinear sde 



In the Section VI we have shown that the herding model 
of return in the financial markets leads to the class of 
stochastic differential equations, whose general form is given 



by SDE ( 18 i. This class of stochastic differential equations 
reproduces power law statistics, namely the probability 
density function and the spectral density, of return and 
trading activity in the financial markets. The burst statistics 
of the financial markets are also very important for the risk 
management and would serve as an additional criteria to 
determine the model consistency. In this section we provide 
some initial results of burst statistics generated by the SDE 

We define a burst as a part of the time series lying above 
the certain threshold, hj. In Figure |7] we present an example 
burst of the simple bursty time series, I{t). Evidently a burst 
as itself can be described by its duration, T = t2 — ti, 
maximum value, Imax, and burst size, which we define as 
an area above the selected threshold yet bellow time series 
curve (highlighted by x pattern in the Fig. |7|i, S. 

There is a well established passage, or alternatively hit- 
ting, time framework, which is frequently used to tackle 
practical problems in both mathematical finance ll34l and 
physics f26l , 1149 II . One can also apply this framework to 
understand the burst durations, T. Interestingly enough we 
can consider the first hitting time of the stochastic process 
starting infinitesimally near the hitting threshold as the burst 
duration itself, T. 

Brownian motion, geometric Brownian motion and Bessel 
process are highly applicable models (for examples of the 
application in the mathematical finance, see Ii34i ) for which 
hitting times statistics are known. The Bessel process, 

^W.. (45, 



di? = 



R 



is one of the most interesting as some prominent math- 
ematical finance models can be transformed to a similar 
form. In order to simplify further handling of the Bessel 
process it is convenient to introduce — ^ — 1, which 
is known as the index of the Bessel process. While N 



is also frequently retained and mentioned as it bears an 
actual physical meaning - the Bessel process is an Euclidean 
norm, length of the vector, of iV-dimensional Brownian 
motion, which starts at the origin. Note that for > 1, or 
alternatively i/ > —0.5, R tends to diverge towards infinity. 

In our case the Bessel process is of high interest as by 
using the Lamberti transform defined as 



i-y^ z{y) 



1 



(77 - 1)2;"- 



(46) 



we can reduce a general class of SDE ([18} to the Bessel 
process. 



dz — { 11 



2 z 



\{'n-t) ■ corresponding dimension of 



with index 1^ 

the Brownian motion is given byA^ = 2(i^ + l) = ^5y- 

Let us assume that a burst starts at time t^, with = 
y(to) slightly exceeding the selected threshold, hy. By 
definition the burst lasts until y(t) crosses hy once again, 
but now from the above. Equivalently, in the terms of Bessel 
process the burst lasts until at a certain time, t, the z process 
crosses the boundary h^^ = £{hy) from the below, while the 
starting position, zq ~ z{to), which in the terms of Bessel 
process is below the threshold, zq = i{yo) < h^- 

Consequently by choosing zq arbitrarily close yet below 
hz, we can obtain an estimate for the burst duration, T, in 



terms of the hitting times of the Bessel process, r 



'T: — ml 

t>to 



[t, z{t) > /i,}, 

Zq < e, 



(48) 



where e is an arbitrary small positive constant. As given in 
II5OJ . the following holds for < zq < 



'v-2 °° 

fc=i 



■^0 



■ exp 



'2hi 



(49) 



where pf^h (0 ^ probability density function of the 
hitting times at level of Bessel process with index v 
starting from zq, is a Bessel function of the first kind of 
the order i', while ^ is a fc-th zero of J^. 

We have to replace p^o'/i (^) density function regarding 
hz to avoid the self-evident convergence of p^^J/^ (t) (for t > 
0) to zero, when zq h^- This is achieved introducing the 
probability density function p^j^^ (t) as a probability density 
function of the burst duration 



Phjit)^ lim — , 

zo-s-fij Hz — Zq 



(50) 
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Figure 6. The broad spectra of Hurst exponents, h{q), obtained from time series obtained by solving j32| . The model parameters were set as follows: 
ei = 1, £2 = 2 — a, « = (red squares), 1 (green circles) and 2 (blue triangles). 




Figure 7. Time series exhibiting bursty behavior, I{t). Here hj is the threshold value, above which bursts are detected, ti is the three visible threshold 
passage events, Imax is the highlighted burst's peak value. Burst duration we define as: T = t2 — ti. 



where we have selected the threshold at level and i/ is 
the original model parameter To evaluate this limit we have 



to expand (^j^ji,,kj neai- 



zq _ 1 . 



■ WJi^Ui^m) ~ ii^MJi+Aii^.k)] = (51) 



By using this expansion we can rewrite (50i as: 



k=l 



2Ki 



(52) 



here Ci is a normalization constant. By taking a note that 
i^ k are almost equally spaced, we can replace the sum by 
integration 



G 



43/2 



(53) 



From the expression above follows that the probability 
density of the burst durations in the time series obtained 



by solving SDE ( [181 ) can be approximated by a power law 
with exponential cut-off. Or mathematically 



t 



-3/2 



for t < 



2ht 



for t > 



2ht 



(54) 
(55) 



This result is in agreement with a general property of one 
dimensional diffusion processes presented in 1*491, namely 
that the asymptotic behavior of first hitting times is a power 
law i^^/^ irrespectively of the nature of stochastic one 
dimensional process or the actual mathematical expressions 
of the Langevin or the Fokker-Plank equations. The expo- 
nential cutoff for longer burst durations can be explained 
by the direction preference of the Bessel processes (note 
the positive drift term in case of > 1, or alternatively 
V > —0.5). The actual empirical data, as shown in Fig. [8] 
(b), also has the predicted asymptotic behavior, though the 
inconsistence in fitting is clearly higher than for the model's 
probability density Fig. |8] (a). 

Our empirical data set includes all trades made on NYSE, 
which were made from January, 2005 to March, 2007 and 
involved 24 different stocks, ABT, ADM, BMY, C, CVX, 
DOW, FNM, GE, GM, HD, IBM, JNJ, JPM, KO, LEY, 
MMM, MO, MOT, MRK, SEE, PEE, T, WMT, XOM. We 




have used one hour window moving average filter on empir- 
ical one minute return series. As we consider the model to 
be universal, i.e. applicable towards the modeling of varying 
financial markets and stocks, we can consider each stocks' 
time series as a separate realization of the same stochastic 
process. Time series are first normalized and later averaged 
over the whole set. We back this approach by recalling that 
in O we have shown that the more sophisticated versions of 
( [24) 1 may be well used to model absolute return of different 
stocks from NYSE and Vilnius Stock Exchange. 

There are numerous reasons for the observed inconsis- 
tence in fitting of empirical data Fig. |8](b). Firstly, we were 
unable to remove intra-day pattern from the time series. 
But the main reason is that the simple stochastic model, 
driven by ( fTS] !, is unable to reproduce the full complexity of 
empirically observed spectral density. In order to reproduce 
the correct, fractured, shape of the spectral density one must 
use double stochastic model, driven by a more sophisticated 
version SDE (24i, |3|. Nevertheless, derived equations for 
the burst duration distribution ( [52| and ( [53] l of the general 
process (fTSjl are in agreement with empirical time series 



of return. This provides one more argument for the further 
development of stochastic models based on herding behavior 



of agents and nonlinear SDE (18 1. 



IX. Conclusions and future work 

Reasoning of stochastic models of complex systems by 
the microscopic interactions of agents is still a challenge 
for researchers. Only very general models such as Kirman's 
herding model in ant colony or Bass diffusion model for 
new product adoption have well established agent based 
versions and can be described by stochastic or ordinary 
differential equations. There are many different attempts 
of microscopic modeling in more sophisticated systems, 
such as financial markets or other social systems, intended 
to reproduce the same empirically defined properties. The 
ambiguity of microscopic description in complex systems 



is an objective obstacle for quantitative modeling. Simple 
enough agent based models with established or expected 
corresponding macroscopic description are indispensable in 
modeling of more sophisticated systems. In this contribu- 
tion we discussed various extensions and applications of 
Kirman's herding model. 

First of all, we modify Kirman's model introducing 
interevent time T{y) or trading activity 1/T{y) as func- 
tions of driving return y. This produces the feedback from 
macroscopic variables on the rate of microscopic processes 
and strong nonlinearity in stochastic differential equations 
responsible for the long range power-law statics of financial 
variables. We do expect further development of this approach 
introducing the mood of chartists as independent agent based 
process. 

Nonlinear SDE's derived from the agent herding model 
generate multifractal time series. This gives more confidence 
in the modeling of multifractal series observed in financial 
markets. We derive PDF of burst duration for the basic 



form of nonlinear SDE (18 1. This is in agreement with 
empirical time series of return. Further investigation of burst 
statistics in financial markets in comparison with analytical 
results from nonlinear SDE is ongoing. This would serve as 
an independent method to adjust model parameters to the 
empirical data. 

One more outcome of Kirman's herding behavior of 
agents is one direction process - Bass diffusion. This simple 
example of correspondence between very well established 
microscopic and macroscopic modeling becomes valuable 
for further description of diffusion in social systems. Models 
presented on the interactive web site [lOJ have to facilitate 
further extensive use of computer modeling in economics, 
business and education. 
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